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ABSTRACT 

A  mathematical  model    for   the   delivery   of   conventional   bombs   by 
stick   bombing    has    been  developed.      This   model   considers    bomb    type, 
number   of   bombs,    pilot   miss-distance   distribution,    and   dimensions 
of  a    rectilinear   target.      It   may  be    used    to  estimate    the   effects    of 
bomb    interval   and   approach   angle    upon   the    single    pass    destruction 
probability    of   a    target. 

The    significant   variables    of    the    stick   bombing    operation  are 
identified,    and    functional   relationships   are   developed   among    these. 
The    model    is    used    to   examine    the    problem  of  determining   an   optimal 
approach  angle   and  bomb   interval    for   stick   bombing   bridges   and 
railroads . 
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I.   STATEMENT  OF  PROBLEM 

During  the  Vietnam  conflict  the  U .  S.  Navy  was  given  the  primary 
responsibility  of  slowing  the  flow  of  supplies  from  North  Vietnam  to 
South  Vietnam.   This  task  was  partially  accomplished  by  cutting  rail- 
roads, roads,  and  bridges  on  all  North-South  avenues  of  travel  in 
North  Vietnam.   The  major  weapons  systems  used  were  light  jet  air- 
craft* carrying  conventional  general  purpose  bombs.* 

Since  this  type  of  combat  bombing  had  last  been  utilized  during 
the  Korean  War  and  in  fact  had  only  been  practiced  as  a  secondary 
mission  (the  primary  mission  of  the  squadrons  involved  had  been 
nuclear  strike)  limited  knowledge  and  know  how  existed  on  how  best 
to  attack  these  targets.   Furthermore  the  available  publications  at 
that  time  failed  to  provide  adequate  answers.   In  many  cases  the 
publications  discussed  targets  of  vastly  different  dimensions  and 
construction  than  those  encountered  in  North  Vietnam  and  the  recom- 
mendations advanced  in  these  publications  failed  to  fully  incorporate 
the  advantages  of  stick  bombing.* 

The  pilots,  therefore,  had  to  develop  their  own  estimate  of  the 
best  tactics  to  use.   These  tactics  were  largely  based  on  "seat  of 
the  pants"  estimates  and  had  little  if  any  analytical  evidence  to 
support  them.   It  is  the  purpose  of  this  thesis  to  provide  the 
necessary  analytical  techniques. 


*A11  terms  with  asterisks  are  defined  in  Appendix  E 


In  order  to  discuss  the  tactics  developed,  the  variables  available 
to  the  pilot  must  first  be  delineated.   Once  the  target  has  been 
selected  the  pilot  has  some  control  over  the  following  variables: 
(1)  the  number  of  airplanes  used  to  attack  the  target,  (2)  the  actual 
order  and  interval  between  airplanes  in  the  attack,  (3)  the  number  of 
passes  per  airplane  over  the  target,  (4)  the  number  and  type  of  weapons 
that  the  airplanes  carry,  (5)  the  approach  angle*  to  the  target,  (6) 
the  interval*  between  weapons  and  hence  the  stick  length*,  (7)  the 
maximum  and  minimum  altitudes  of  releasing  the  weapons,  (8)  the  dive 
angle*  and  speed  in  the  dive. 

Although  tactics  varied  somewhat  from  airwing*  to  airwing  and 
even  squadron*  to  squadron,  an  accurate  description  of  the  generally 
accepted  typical  tactic  used  to  cut  either  a  road,  railroad,  or  bridge 
is  as  follows :   for  a  given  target  3  or  4  airplanes  loaded  with  4  to  6 
500  lb.  general  purpose  bombs*  were  sent  as  a  flight*.   The  aircraft 
proceeded  in  formation  to  the  target  vicinity  whereupon  the  flight 
leader  identified  the  target  and  rolled  in*,  followed  in  approximately 
5  seconds  and  with  a  different  approach  angle  by  the  second  airplane 
and  then  in  a  similar  manner  by  the  third  airplane,  etc.   Each  pilot 
made  an  individual  tracking*  dive  attempting  to  dive  at  a  predeter- 
mined dive  angle  and  speed.   At  a  preselected  altitude  each  pilot 
released  all  of  his  bombs  in  a  stick  of  fixed  and  equal  bomb  inter- 
vals and  then  departed  the  target  area. 


Joint  Munitions  Effectiveness  Manual  (Air  to  Surface)  Weapons 
Effectiveness,  Selection,  and  Requirements,  1A-4,  Unclassified, 
5  April  1968. 
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The    questions    immediately    raised    when   examining    these    tactics 
are:      for   a   given  bomb    load,   what    are    the   effects    of   changing    the 
approach   angle    on   the    probability   of  a    single    pass    hit*   (Psph)? 
What    does    the    interval   between   bombs    depend    on   and    how   can    it    be 
set    so   as    to  maximize    the    Psph?     What   effects    do   target   dimensions 
have    on    the    tactics    developed?     And    lastly   since    pilot    skill   varies 
from   pilot    to   pilot,    how  does    changing    the    pilot   miss-distance 
distribution*  affect    the    overall    results? 

In   attempting    to  answer   these   questions,    several    of    the    pilot 
controlled    variables   are   assumed    fixed    throughout    the    investigation. 
In   the    first    place    it    is    assumed    that    the    airplanes   attack    separately 
and   make    individual    tracking    runs;    the   dive    angle,    release    speed, 
maximum   and   minimum   release    altitudes    are    invariant    for  all    tactics 
investigated;   and    lastly  all  weapons   are   dropped   on   a    single    pass. 
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II.      DEVELOPMENT   OF   THE   GENERAL  MODEL 

To  answer    the   above   questions    relating    to   tactics,    an   operational 
model    of    the   weapon   system  will   be    developed    in   the    following    se- 
quential  manner:      analyze    the    operation  and    then   translate    these 
results    into  a   mathematical  model. 

In   analyzing    the    operation  of   the  weapon  system   (light    jet   attack 
aircraft    carrying   conventional  general   purpose   bombs)    it   seems 
appropriate    to  address    the    following    topics  : 

a.  What   are    the    significant    parameters    of  weapon  system 
performance    in   its    operational  environment? 

b.  How   is    system   performance    related    to   these    parameters 
(what    functional   relationships)? 

c.  Parameters    used    in  any  mathematical  model   subsequently 
developed   should   be    capable    of  being   estimated  with 
empirical  data. 

The   above    operations   analysis   will   be    used    to  generate   a 
mathematical   model.      This    model  will   be    used    to  gain    insight    into 
the    performance    of   the    system  and   determine    how   to    optimize    system 
effectiveness   according   to   the   criterion   of   single    pass   destruction 
probability. 

Thus,   we   must    start   with  an  analysis    of   the   bombing    operation. 
The    entire    sequence    of    operations    from   target    detection    to  munition 
detonation  was    considered.      System  performance   was    partitioned    into 
the    following    subsystems : 
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a .  Acquisition 

b.  Delivery 

c.  Effects 

Acquisition  was    not    considered    further    since     it    is    a    factor    held 
constant    in    this    investigation.      System  delivery   capability  was    further 
subdivided   as    follows : 

a.  from  detection   to   bomb    release, 

b.  from   bomb   release    to  detonation. 

The    first   aspect    is   discussed    in   the    section,    "Analysis    of  Aircraft 
Attack    from   Roll-in   to  Weapon   Release",   while    the    second    is    discussed 
in    the    section,    "Analysis    of   Pilot  Miss-distance   Distribution".      These 
various   aspects    of   the    bombing  attack  will    now  be   discussed.      It   should 
be    noted    that    the   analysis   was    based    not    only  upon    the   author's    3-1/2 
years    in  a    light    jet   attack    squadron  during  which    time   many  missions 
were    flown  against    targets    of   this    type   but   also   upon  discussions   with 
many   pilots    of   similar   backgrounds. 

A.      ANALYSIS   OF   PILOT  MISS -DISTANCE   DISTRIBUTION 

In   the    first    place,   miss-distance   was    defined   as    the    error   between 
the   actual   point    of    impact    of   the   bomb   center   and    the    intended    point 
of    impact    of    the   bomb   center.      An   orthogonal   coordinate   system  was 
defined   along   and    perpendicular   to    the    flight    path   of   the   aircraft 
with    its    origin  at    the    intended    point    of    impact    of    the    bomb    center 
or   stick   center    if   more    than   one   weapon  were  dropped .      The   error  was 
then  decomposed    into   two   components:      range   error    (R)    along   the 
flight   path   and   deflection  error    (D)    perpendicular   to    the    flight    path. 
It   was    assumed    that    range   error  and   deflection  error  were    statistically 
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independent   with   a    joint    density    function   of    the    bi-variate    normal 
form , 


fRD    (R'D)    '2^5 


exp 


R  D 


R   - 


hr 


^)  .^)]| 


For   a    random   variable    such   as    range    error  which    is    composed    of   many 
factors    (wind,    target   acquisition,    bomb   release,    etc.),    a    postulated 
normal   distribution   is    reasonable    in    light    of   the   Central   Limit 
Theorem.      This    assumption  was   also   consistent  with  available   empirical 
data  . 

Since    all   pilots   keep   extensive    records    of   bomb   hits   made    during 
practice   bomb    runs    on   instrumented    targets,    pilots   were   assumed    to 
know   their   expected    range   error   and   deflection   error  and    to   be   able 
to   correct    for    it,    that    is   \j.     —  (wL    =  0.      This   assumption   is    obviously 
optimistic   but    it   guarantees    that   all    results    obtained    for   Psph  are 
upper  bounds    and    it   certainly   is    the    ideal    toward  which   pilots   work. 

Defined    in  this   manner,    the    pilot   miss-distance   distribution   can 

be   characterized   by    two   parameters,    namely  aD  and   a    ,    and   by  varying 

R  D 

aD   and   a    ,    the   effect    of   the    pilot  miss-distance   distribution   on 
R  D 

Psph   can  be    observed. 

B.      ANALYSIS    OF  AIRCRAFT  ATTACK  FROM  ROLL-IN  TO  WEAPON  RELEASE 

The    parameters    that   affect    the   attack    from   roll-in   to  weapon 
re  lease    are  : 


"Anderson,   T.    W. ,    An   Introduction   to  Multivariate   Statistical  Analysis, 
p.     11,    Wiley,    1958. 
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Dive  Angle 

Release  Speed 

Maximum  Release  Altitude 

Minimum  Release  Altitude 

Approach  Angle 

Bomb  Interval 

In   order    to   examine   different    tactics    available,    a    standard   bomb 
run   was    defined    and    used    throughout    the    investigation.      That    is    to    say 
dive    angle,    release    speed,    maximum  and   minimum    release    altitudes   were 
all    held    constant.      The    purpose    of  defining    this    standard   bomb    run 
was    to  maximize    Psph    for    this    standard.      Furthermore,    these    parameters 
actually   control    the    pilot   miss-distance    distribution.      Hence,    if    the 
examination    of   different    types    of   bomb    runs    is    desired,    only    the 
relationship  between   these    parameters   and    the    pilot   miss-distance 
distribution   need   be   determined   before    the   analysis    is    performed. 

Approach   angle   was    defined   as    the   acute   angle   between   the    center 
line    of   the    target   and    the   horizontal   projection   of    the    flight    path. 
During    the    investigation   the   approach   angle   was    allowed    to  vary   from 
0    -    90    .      It   was    assumed    that    the   approach   angle   was    selected    prior 
to    the    bombing    run  and    that    the    pilot   attempted    to   make    his    attack 
using    this   approach  angle.      In   order    to  account    for    pilot    error    in 
this    parameter,    a   known   probability  density    function,    g    (9) ,    of 
approach   angle    error    (AAE)   was    assumed. 

Lastly   the    only   constraint    placed   on   the    selection   of  bomb 
interval    (INT)   was    that    for   any  bomb   stick    the    intervals    must   be 
equa 1 . 
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With  these  assumptions,  it  is  possible  to  mathematically  describe 
the  attack  from  roll-in  to  weapon  release  by  three  parameters:   the 
number  of  bombs  to  be  dropped  (n) ,  the  preselected  approach  angle 
(0  ) ,  and  the  bomb  interval  (INT). 

C.   ANALYSIS  OF  MUNITION  EFFECTIVENESS 

The  weapons  and  the  weapon  release  mechanisms  were  assigned  a 
reliability  factor  of  1.   For  a  comparative  evaluation  of  tactics 
this  is  reasonable,  however,  an  analysis  of  the  effect  of  relaxing 
this  assumption  was  made  and  a  discussion  of  it  can  be  found  in 
Chapter  IV. 

A  convenient  method  of  describing  weapons  effects  is  to  define 
a  lethality  function,  L(d),  where  the  lethality  function  is  a 
description  of  the  probable  damage  caused  as  a  function  of  miss- 
distance.   For  example,  for  a  one  dimensional  target  the  lethality 
function  might  take  the  form  of  Figure  1. 


Kill  probability 


i 


Target 


Typical  one  dimensional  lethality  function 
FIGURE  1 
For  this  model  a  "cookie  cutter"  type  lethality  function  was 
assumed.   Although  it  has  been  noted  the  lethality  of  a  fragmenting 
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munition    obeys    an  exponential   decay   with    radial   miss-distance    ,    the 
exploratory    nature    of    this    investigation   seemed    to    justify    the    simpler 
model.      Thus    the    probability   of  killing*   the    target    is    1    if    the   miss- 
distance    is    less    than   or    equal    to  R   and    0    otherwise.      See   Figure    2    for 
the    one   dimensional   case. 


Probable   damage    factor 


One   dimensional    "cookie    cutter"    lethality    function 

FIGURE   2 
By    using    this    type    of    lethality    function    the   weapon   effectiveness 
can  be    characterized   by   one    parameter.      In    this    case    the    parameter    is 
defined    as    the    effective    miss-distance    (EMD)    and    is    the    maximum 
distance    from  any   target    dimension   that    the    center   of    impact    of  a 
weapon  may  be   and    still   produce   a   kill. 

D.   MATHEMATICAL  FORMULATION 

The  problem  has  now  been  reduced  to  investigation  of  the  effect 

on  Psph  of  6  parameters,  namely 

n  =  number  of  bombs  dropped , 

aD  =  standard  deviation  in  range  error, 
R 

O     =  standard  deviation  in  deflection  error, 
INT  =  bomb  interval, 


Ballistics  Research  Lab  Report  697,  Justification  of  an  Exponential 
Fall  Off  Law  for  Numbers  of  Effective  Fragments,  by  H.  W.  Weiss, 
February,  1949. 
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9   =  preselected  approach  angle, 

EMD  =  effective  miss -distance  , 
and  one  arbitrarily  defined  probability  density  function,  the  approach 
angle  error  density,  g  (9)  . 

In  the  first  place,  the  number  of  weapons  (n)  was  assumed  known 
and  given  for  any  attack.   Hence  for  a  given  n 


T    -r 


max  Psph  =    J  J  L(R,D)  fR  D<R,D)  gQ(9)  dRdDd9, 

all  space   all  9 

where  L(R,D)  is  the  lethality  function  set  equal  to  1  when  (R,D)gA*(9) 
(0  elsewhere)  and  A*(0)  is  the  area,  in  the  R,D  plane,  depending  upon 
the  approach  angle,  9,  over  which  the  stick  center  may  range  and  still 
achieve  a  target  kill.   Obviously  A*(9)  is  also  a  function  of  EMD,  INT, 
and  target  dimensions.   In  order  to  describe  A*(9),  however,  the  opti- 
mum bomb  interval  must  first  be  determined. 

As  seen  in  Appendix  A,  by  examining  the  geometry  of  the  situation 
and  the  characteristics  of  the  normal  density  function,  the  optimal 
bomb  interval  can  be  developed  and  is  a  function  of  EMD,  9  and  target 
dimensions.   It  is  shown  in  Appendix  A  that  the  optimal  bomb  interval, 

INT    ,  is  given  by: 
opt 

INT    =  [2  EMD  +  W]/sin  9. 
opt 

In  order  to  perform  the  integration,  it  is  now  necessary  to 

determine  the  boundaries  of  A*  in  terms  of  the  optimal  bomb  interval 

or  in  other  words  in  terms  of  EMD,  target  dimensions  and  9  .   This 

can  be  done  but  due  to  the  complexity  of  the  development  and  the 

differences  in  target  shape,  the  details  appear  in  Appendices  B  and 

C.   Once  the  boundaries  have  been  determined,  however,  and  represented 
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MAX    Psph    =         j   I  J        gg(9)    fR   D(R,D)    dRdDdQ 


in  the  R-D  coordinate  system 

[j 

R,DeA*  All  0 
may  be  numerically  integrated  by  the  method  discussed  in  Appendix  D. 
Thus,  optimal  tactics  may  be  determined  by  varying  the  parameters  on 
which  Psph  depends. 
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III.   APPLICATION  OF  MODEL  IN  THE  DETERMINATION  OF  OPTIMAL  TACTICS 

A.   RAILROAD/ROAD  TARGET 

1 .   Form  of  the  Model 

In  this  case  the  length  of  the  target  is  assumed  to  be  so 

great  that  it  is  effectively  impossible  to  miss  the  target  along  its 

length  dimension.   This  is  equivalent  to  assuming  that  the  railroad/ 

road  is  an  infinitely  long  rectilinear  target.   A  straight  segment  of 

roadway  approximately  18001  long  is  an  example  of  a  target  which  may 

be  considered  to  be  infinitely  long.   In  order  to  show  this  suppose 

that  the  approach  angle  (9)  is  0  and  that  only  one  weapon  is  to  be 

dropped.   Since  the  range  error  is  always  assumed  to  be  larger  than 

deflection  error,  this  approach  angle  causes  the  length  of  the  target 

to  be  most  critical.   Because  99.64$  of  the  probability  of  a  normal 

distribution  lies  between  +  3a,  6a_  will  include  essentially  all  the 

R 

possible   errors    in  the   range   direction.      The    largest   value    of   a     used 

R 

in   this    investigation  was    300    feet;    hence   any   target    1800    feet    or 

longer  will    produce    results    that    do    not   vary   significantly    from 

results    obtained  when  a    target    of    infinite    length    is    used. 

Suppose    that    the    target   has  width    (RRW) ,    length    (L >  ») , 

that    the    distribution   is   bi-variate   normal  with  known  o     and   a_   and 

R  D 


UR   =  ^D   =   °' , 

R  D  aR         aD2 

Define    two   coordinate    systems,    a   X-Y   coordinate   axis    parallel   and 
perpendicular    to    the    target    dimensions    and   a   R-D   coordinate    system 
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(used  to  measure  pilot  error)  parallel  and  perpendicular  to  the 
flight  path.  D         Y 


Coordinate  axes  in  railroad/road  problem 
FIGURE  3 
In  order  to  reduce  the  notation  assume  that  n  =  2  (i.e.,  the  stick 
consists  of  two  weapons).   Furthermore,  for  any  given  approach  angle 
(9),  set  the  bomb  interval  (INT)  for  these  weapons  at  the  optimal 
value 

INT  =  (2  EMD  +  RRW)/sin  9. 
Now,  there  is  a  target  kill  if  the  CS  falls  within  certain  limits. 
See  Figure  4. 


Limits  of  position  of  stick  center  needed 
produce  a  target  kill 

FIGURE  4 
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where  y   =  upper  limit  for  CS  and 
y   =  lower  limit  for  CS . 

These  limits  of  integration  can  be  computed  and  expressed  in 
the  R-D  coordinate  system.  See  Appendix  B,  hence  A*  is  well  defined 
and 

r 

MAX  Psph  =  f  ^(f^D)  gQ(9)  dRdDdO 

R,DeA*  All  9 

can  be  numerically  evaluated  for  each  9  ,  n,  EMD,  target  size,  and 

V  V 

2 .   Sample  Results 

The  results  presented  in  this  section  depict  only  several  of 
the  many  possible  combinations  of  parameters  available.   The  results 
are  presented  in  order  to  illustrate  the  possible  functional  depend- 
ence of  angle  and  bomb  interval.   For  this  computer  run  the  following 
values  of  the  parameters  were  used:   (1)  bomb  load  =  6  (2)  EMD  =  12.5' 

(3)  railroad  width  =  5'  (4)  100'  <  aD  <  3001  (5)  50'  •:  a  <:  150'. 

R  D 

Since    no  data   was    available    to  estimate    the   probability  density 
function   of  approach  angle   error,    the   author  assumed   a   discrete   density 
function   that,    based    on  his    experience,    seemed    representative.      It  was 
defined   as    follows    for   every   preselected  approach  angle    9 


gQ(9)=.l        00   -    9 


0 
o 


-.1  80-6° 

=    .2  o„   -    3° 

=    .2  9     +  3° 

-  .1  %  +  6° 
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,.1        %+9° 


where    9      =   pre-planned   approach   angle. 

The    optimum  approach  angle   was    found    to   be    8.62      but    Psph    for 
this    angle    did    not    vary   appreciably    from    the   value    obtained   when 
9   =    10    .      Notice    that   at    10    ,    the    rate    of   change    of    the    bomb    inter- 
val   is    relatively    large;    hence,    the   assumption  about    proper    function- 
ing   of    the    bomb    release    equipment    might    be    critical.       On    the    other    hand, 
when   this    particular   distribution    for  Approach  Angle   Error    (AAE)   was 
assumed,    the    results    obtained   at    9      =    10     are    not    too  different    from 
those   when   no   error  was    assumed.      The    same    type    of    results    could   have 
been   obtained   by  assuming    that    there   were    no  approach   angle   errors 
but    that    stick    length  was   governed   by  an  error   distribution  since   both 
types    of   errors    cause   essentially    the    same    results    (i.e.    a   bomb    inter- 
val   is   being   used    that    is    not    optimal    for    the    particular   approach 
angle).      Furthermore,    notice    that   by   assuming    that   stick    length   and 
approach  angle   are    free    from  error   although    not    necessarily   causing 
erroneous    results    tend    to    inflate    the   values    of   Psph   and    to   cause    the 
maximum  Psph    to   occur   at    too   shallow  an  angle. 

It    is   most    important    to   remember    that    the   values    of   Psph  were 
all   determined   by   calculating   and    using    the    optimum  bomb    interval    for 
the   approach   angle    or    the   E    (Approach  Angle)    in   the    case    of  AAE   and 
these    values   would    be    considerably   different    if    for   example^  the    bomb 
interval    for   9   =   80     were    used    in  calculating    the    Psph    for   9   =    10    . 

Notice    that    Psph   is   almost    constant    for   any  given   a     through- 
out   the    investigated    range    of   a    .      It    is    only  when  a     =   50'    that 

R  D 

appreciable  fall  off  occurs  as  aD  increases.   This  indicates  that  if 

R 

the  standard  deviation  of  pilot  error  in  deflection  is  greater  than 
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50 ' ,    effort    should   be    expended    to   reduce    it   prior   to  attempting    to 

reduce    a    .      For   example,    suppose    G     =    100'    and    a     =  300',    then   if 
R  UK 


O     is    reduced   by   50^,    a     held   constant,    Psph   increases    from    .60   to 
D  R 

.83,    but    if  a      is    reduced   by   50$,    <j     held   constant,    Psph   only   increases 
R  D 

from    .60    to    .62. 

Also  it  should  be  noted  that  approaching  the  target  exactly 
along  the  center  line  of  the  target  produces  the  lowest  Psph.   This 
of  course  is  caused  by  the  fact  that  at  0  =  0  all  advantages  of 
stick  bombing  are  lost.   Finally  when  0  is  larger  than  45  ,  two 
interesting  facts  appear:   (1)  the  bomb  interval  is  almost  constant 
and  (2)  when  EAA  is  assumed  to  be  present,  the  results  compare  almost 
exactly  to  those  obtained  when  EAA  is  assumed  absent.   These  results 
can  be  interpreted  to  mean  that  if  errors  in  approach  angle/bomb 
interval  are  thought  to  be  large  and  unpredictable  (i.e.,  the  error 
distribution  cannot  be  estimated)  perhaps  more  consistent  results  can 
be  obtained  by  using  approach  angles  greater  than  45  . 

B.   BRIDGE  TARGETS 

1.   Form  of  the  Model 

Because  of  its  finite  length,  the  area  over  which  a  kill 
occurs  is  somewhat  different  than  the  infinite  strip  involved  in 
the  railroad  problem.   In  order  to  examine  this  area,  define  an 
X-Y  coordinate  system,  parallel  and  perpendicular  to  the  target 
dimensions . 
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GRAPH   II 
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GRAPH  III:   Optimal  Bomb  Interval  vs.  Approach  Angle 
EMD  ■  12.5' 
n  -  6     Railroad/Road  Target  5'  Wide 
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^x 


Bridge  target 
FIGURE  5 
Now  increase  the  target  dimensions  by  a  quantity  representing  the 
effective  miss-distance  for  the  weapon  in  question.   Call  the  new 
dimensions  effective  bridge  length  (EBL)  and  effective  bridge 
width  (EBW). 


Effective  bridge  dimensions 
FIGURE  6 
Notice  that  the  shaded  areas  are  approximations  to  the  actual  dis- 
tance away  from  the  target,  that  the  centers  of  impact  may  be  and  still 
produce  a  kill.   The  small  increase  in  accuracy  obtained  by  requiring 
the  boundaries  to  be  exact  is,  however,  far  outweighed  by  the  com- 
plexity introduced  in  the  mathematics  involved;  hence,  the  approxi- 
mation is  used  throughout  the  development. 
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By   defining    the    target    in   this   manner,    a    target   kill    is 
produced   whenever    the    center    of    impact    of   a   weapon    falls   within    the 
rectangle    defined    by   EBL  and    EBW. 

The    first    problem    is    to    recognize    that    there    are    two   basically 
different    approach  angles   available    in  attacking    the    target:      (1) 
less    than    the    critical   angle,    and    (2)    equal    to    or    greater    than    the 
critical   angle   where    the    critical   angle    is    defined   as: 


Critical   angle    =  p   =   tan 
Y 


-1 


EBW 
LEBL-T. 


Critical  approach  angle 
FIGURE  7 
Both  types  of  attack  are  developed  but  since  they  differ  only  in  the 
description  of  the  integration  limits  only  the  case  when  9  >  @  is 
discussed  below.   For  ease  in  notation  assume  n  =  2.   The  optimal 
bomb  interval  is  determined  in  exactly  the  same  manner  as  for  the 
railroad/road  problem.   That  is 

INT  =  (2  EMD  +  BW)/sin  9  =  EBW/  sin  9. 
Next,  the  question  of  where  the  center  of  the  stick  can  fall  and 
still  produce  a  kill  must  be  answered.   Once  this  has  been  deter- 
mined and  transformed  into  the  R-D  coordinate  system,  see  Appendix 
C,  the  problem  is  essentially  the  same,  that  is 
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MAX 


Psph    =       //        fR  D(R'D)    dRdD._ 


proper 
polygon 

Since    the   distribution   of  approach  angle   error   produced   only 

minor   variations    in   the    results    in   the    railroad    problem,    in   this    case 

g    (0)   was   assumed    to   have   all   of   its    probability  mass   at   9    ,    hence    its 
y  u 

effect   does    not   appear   in   the    integral    for   Psph. 

By   using   the   numerical   method   of  Appendix  D   to  evaluate    the 

above    integral,    the   maximum  Psph    for  any  target   dimensions,    a    ,    a    , 

K        D 

n,    9  can  be   determined. 
2 .      Sample   Results 

Results   are    once   again   presented    for   only   one    representative 
combination   of   parameters   and  all   comments   made   apply   only   to   this 
set    of   parameters.      The    parameters    used  are    (1)   n  =  6   bombs,    (2) 
EMD  =  5',    (3)    40'    <  bridge    length      <     200',    (4)   bridge   width   = 
10'    and  20',    (5)    100'   •:  a,  <  300',    (6)   a     =  50',    100'. 

The   graphs   are   rather  self-explanatory  but    perhaps    several 
facts   are  worthy   of  comment.      In   the    first    place   until  a   bridge 
length   of  approximately    100'    is    reached,    Psph    is    relatively   unaffected 
by  approach  angle.      Remember,    however,    that    (1)   each   Psph   is    computed 
using    the    optimum  bomb    interval   for    the   approach  angle    in  question, 
and    (2)    selecting  a   bomb    interval   that    is    not   compatible  with   the 
approach   angle  will   produce    lower   results.      Once   again   stick    length 
and   bomb   interval   change    rapidly  at    the   shallow  angles   and   are 
relatively   constant    for   9  >  45    .      These    two   facts    point    toward 
approach  angles    in  excess    of  45      if    large   approach  angle   error   or 
bomb    interval   error    is   anticipated. 
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Notice  also  that  as  the  bridge  dimensions  increase  the 

difference  between  the  Psph  for  a  pilot  miss  distribution  of 

G  =  150',  an  =   50'  and  a  pilot  miss  distribution  of  o     =  300', 
R  D  R 

0_  =  100'  increases.   For  example  when  the  bridge  length  is  increased 

from  100'  to  200'  Psph,  for  pilot  gd  =  150',  a^  =  50',  increases  66$, 

R  D 

while  Psph  for  a  pilot  o     =   300',  o_  =  100'  increases  only  59$. 

R  D 
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GRPAH  IV:   Maximum  Probability  Single  Pass  Hit  vs. 
Approach  Angle 

EMD  =  5'  n  =  6 

Bridge   Target  Width,    10',    20';   40'    <    length  <  200' 
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GRAPH  V:   Maximum  Probability  of  Single  Pass  Hit  vs. 
Approach  Angle 

EMD  ■  5*  n  =  6 

Bridge  Width  10',  20';  40'  <  Length  <  200' 
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GRAPH  VI:   Maximum  Probability  Single  Pass  Hit  vs. 
Standard  Deviation  in  Range 
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GRAPH   VII 
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GRAPH  VIII:   Optimal  Bomb  Interval  vs.  Approach  Angle 
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IV.      CONCLUSIONS 

In   conclusion    it    is    evident    that    for   any  given    target    dimensions, 
bomb   size,    number    of   bombs,    and    pilot    miss-distance    distribution, 
there    is    definitely   an  approach   angle    and    bomb    interval    that    maximizes 
the    probability    of   a    single    pass    hit. 

Furthermore    it    has    been   shown   that    for  all    targets    considered: 

(1)  the  optimal  bomb  interval  is  independent  of  the  pilot 
miss-distance  distribution  and  can  be  computed  by  the 
formula  : 

INT   =    [2    EMD  +  W]/sin   9 

where        EMD   -  effective   miss-distance    of    the   weapon 

W     =   target  width 

9      =   preselected   approach   angle, 

(2)  optimal   stick    length 

SL  =  (n-1)  [2  EMD  +  W]/sin  9 
where    n  =  number  of  bombs  dropped, 
and 

(3)  although,  in  general,  no  closed  form  exists  for  the 
optimal  approach  angle,  for  all  cases  investigated 
the  optimal  angle  was  in  the  neighborhood  of  10°-30°. 

This  latter  angle,  however,  is  critical,  since  Psph  is  quite  sensitive 

to  it.   For  example  in  the  railroad  target  problem,  while  approaching 

the  target  from  approximately  10   produced  a  maximum  Psph,  approaching 

from  0   produced  a  minimum.   Likewise  for  the  bridge  target  problem, 

approaching  the  target  from  an  angle  slightly  greater  than  the  diagonal 

produced  a  maximum  Psph,  but  when  the  approach  angle  became  less  than 

the  diagonal,  the  Psph  was  greatly  reduced. 
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It  was   also   shown   that    the    rate    of   change   of    the    optimal   bomb 
interval  was   much    larger    for   small  approach  angles    than   for  approach 
angles   greater    than  45    .      In   fact    for  approach  angles   greater   than 
45    ,    the    optimal   bomb    interval  was    almost    constant.      These    facts 
indicate    that    it  may  be   better   to   use   a    non-optimal  approach  angle 
greater   than   45      if    large   approach  angle/bomb   interval   errors   are 
anticipated . 

For   the   railroad/road    problem,    it  was    shown  that    the    optimal 

approach  angle   also  depends    on   the  maximum  and  minimum  release 

altitudes   and   the   dive   angle    (a)   by   the   relationship 

ir(ALT          -  ALT    .    )   cot    (a)-, 
-  IT  max mm 

opt   5S  L(n-l)    (2   EMD  +  RRW)  J 

where   a  -  dive   angle, 

RRW  =  railroad  width. 
Once   again,    however,    approach  angle    is    independent    of   pilot  miss- 
distance   distribution.      Furthermore    for    the   cases    investigated, 

Psph  was    relatively   constant    for  all  values    of  a_,   when   a     was 

R         D 

greater  than  50'.   This  indicates  that  effort  should  first  be 
expended  to  reduce  the  error  in  the  deflection  direction  if  o_  is 
greater  than  501.   Lastly,  it  was  shown  that  if  approach  angle  error 
was  introduced,  only  minor  variations  of  the  results  were  produced, 
and  for  9  >  45  ,  this  variation  was  negligible. 

In  the  bridge  problem,  it  was  shown  that  for  bridges  less  than 
100  feet  in  length,  Psph  is  not  sensitive  to  approach  angle  as  long 
as  the  optimal  bomb  interval  is  used  for  the  preselected  approach 
angle.   As  the  bridge  length  increased  the  results  asymptotically 
approached  those  of  the  railroad  problem.   It  should  be  noted, 
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however,    that    it   was    impossible    to   calculate    the    optimal   approach 
angle    prior    to    performing    the    numerical    integration,    but    in  all   cases 
the    optimal   approach  angle  was    slightly  greater    than   the    angle    formed 
by    the    diagonal. 

One    last   comment   must    be   made    about    the    results    derived    from   the 
models   and    the    conclusions   drawn   therefrom.      It   has    been  assumed    that 
a    pilot   will   use    the    optimal   bomb    interval    for    the    preselected   approach 
angle    or  expected   approach  angle.      Either   the   derived    results   and/or 
the   conclusions   drawn   therefrom  may   cease    to   be   applicable    if    the   bomb 
interval   and   approach  angle   do   not    satisfy   the   equation: 

INT  =    [2   EMD  +  W]/sin   0    . 
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APPENDIX  A 

OPTIMAL  BOMB  INTERVAL 

Before  the  formulation  of  an  optimal  bomb  interval  can  be 

developed,  some  discussion  as  to  how  stick  bombing  affects  pilot 

error  is  in  order.   Stick  bombing  is  based  upon  the  fact  that 

during  the  tracking  dive,  the  pilot  is  able  to  release  his  weapons 

at  programmed  intervals.   These  weapons,  at  least  theoretically, 

fall  in  a  "straight  line"  (as  modified  by  gravity)  and,  except  for 

wind  and  ballistic  effects,  directly  along  the  flight  path.   Implicit 

in  this  assumption  is  the  fact  that  any  wind  or  ballistic  error  affect 

all  bombs  uniformly  and  thus  do  not  disturb  the  "straight  line"  effect 

of  the  tactic.   Furthermore,  since  during  the  tracking  run  the  pilot 

attempts  to  maintain  an  essentially  straight  line  dive  at  least 

during  the  release  of  the  weapons,  any  deflection  error  will  affect 

all  weapons  uniformly.   It  is  evident  then  that  stick  bombing  does 

not  affect  deflection  error  but  does  have  an  affect  on  range  error. 

In  other  words,  by  changing  the  stick  length  (i.e.,  the  bomb  interval), 

for  any  given  pilot  miss -distance  distribution  the  probability  of  a 

single  pass  hit  (Psph)  changes,  and  this  change  depends  upon  the 

relationship  between  aD,  stick  length  (SL) ,  bomb  interval  (INT),  and 

R 

of  course  target  size  and  bomb  effectiveness  (EMD). 

Although  bomb  interval  and  a     are  needed  in  order  to  calculate 

R 

Psph,  it  is  possible  to  show  that  development  of  the  optimal  bomb 

interval  (criteria:   maximize  Psph)  does  not  in  fact  depend  on  c 

R 
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but  only  on  the  target  dimensions,  bomb  type  (i.e.,  Effective  Miss 

Distance,  EMD)  and  approach  angle. 

From  the  definitions  of  Stick  Length  (SL)  and  Bomb  Interval  (INT), 

it  is  clear  that 

SL  =  INT  (n-1) 

where  n  =  number  of  bombs  dropped.   Now  suppose  that  the  target 

(either  railroad,  road,  or  bridge)  has  width  (W)  and  length  (L) ,  and 

that  the  EMD  of  each  weapon  is  M  feet.   See  Figure  8. 
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Target  and  EMD  of  bomb 
FIGURE  8 
Assume  furthermore  that  the  approach  angle  (0)  is  known.   This 
assumption  can  be  viewed  in  two  different  ways  :   (1)  it  applies  if 
the  pilot  preselects  his  approach  angle  and  makes  his  approach  at 
that  angle  without  error  or,  (2)  it  applies  if  the  pilot  preselects 
his  approach  angle  and  attempts  to  approach  the  target  at  that  angle 
but  has  an  error  of  approach  angle  (EAA)  . 

Now  let  us  assume  that  the  number  of  bombs  in  the  stick  is  2, 
(n  =  2).   Once  the  development  for  two  bombs  has  been  completed  the 
results  can  readily  be  extended  to  sticks  of  2 ,  3,  ...,  n  bombs. 
See  Figure  9. 
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Angle  of  approach  and  stick 
FIGURE  9 
Define   CS  =  geometric  center  of  the  stick  with  coordinates 

v  CS'   CS' 
d  =  distance  from  CS  to  Center  of  Impact  (CI)  of  the 

first  weapons  on  either  side  of  CS 

Clearly  2d  =  INT 

and  (n-l)(INT)  =  SL  =  INT. 

Now  suppose    that  CS    falls    upon   the    center    line   of   the    target.      See 

Figure  10. 


>• 


X 


Bomb  stick  superimposed  on  target 
FIGURE  10 
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Define    a   X-Y  coordinate    system  along   and    perpendicular    to    the    center 
line    of    the    target.      Then  a    necessary   condition    for   a   kill    is    to 
require    the   CI    of   either  weapon   to   be   within  +    [EMD  +  W/2]    (as 
measured    in   the    Y-direction)    of   the    target    center    line.      See   Figure 
11. 


^X 


CI  limits  needed  to  produce  a  target  kill 
FIGURE  11 
It  is  not  a  sufficient  condition  because  the  error  in  the  X-direction 
may  be  so  large  as  to  place  the  entire  stick  off  the  target  area. 
This  fact,  however,  does  not  affect  the  development  of  the  optimal 
stick  length  because  the  expected  aim  point  is  the  geometric  center 
of  the  target  area. 

Define:   b  *  M  +  W/2, 
then  the  necessary  condition  for  a  kill  is 

"b  <  fYcs  +  dl  sin  e  £  b 
or 

-b  <  [Y    -  d]  sin  0  <  b. 

The  first  inequality  produces 

b 


I  '  d  -  YCS  -  sin  9 
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and  the  second  inequality  produces 

-   -b  Q  +  d  <  Y_c  <   .b  r  +  d. 
sin  0     —  CS  —  sin  9 

There    are    two    possible   conditions    to    investigate,   — ; -<  d,    and 

sin  9 


sin  9 


>  d 


In  the  first  place,  suppose  — : <   d.   See  Figure  12.   For  this 

rr  sin  9 

case    two   bombs    never   simultaneously  hit    the   target   and    the    probability 
of  a    single    pass   hit    (Psph)    in   the   range   direction   is   given  by 


-   d  + 


sin  9 


Psph   = 


fR   (R)   dR  + 


d  + 

r 


sin   9 


J 


fR   (R)   dR. 


sin  9 


sin  9 


X 


Stick,    b/sin   9  <  d 

FIGURE    12 

Since    f      (R)  ~  NORMAL   (0,    a_),    normalizing   this   expression  produces 
K  R 


Psph   „  |  Q.   d+b/slng^  _   $  ^_   d-b/sin9>)  +  $  ^d+b/sinC^  .  (|  Q 


d-b/sin9^ 
JR 


which  because  of  the  symmetry  of  the  normal  distribution  reduces  to 
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Psph   =  2 


" $  (  d+b/sin9^  _    .  /  d-b/sinQN  "1 


R       '  "  R 

Now  suppose   d    increases,    then   Psph  decreases    because 

when  a '    >  a . 
See   Figure    13 


When 


a'+c 


Normal  density  function 
FIGURE  13 

<  d,  the  Psph  decreases  as  d  increases 


sin  9 
Now  suppose  b/sin  9  >  d.   See  Figure  14. 


Stick,  b/sin  9  d 
FIGURE  14 
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Now 

d  +  b/sin   9 

Psph   =      I  fR(R)    dR 

-   d    -   b/sin   0 
and   by    normalizing 

Psph  =  2  [i  (d  +  b(°lp9)  -   *    (0)]. 

R 

Clearly  as  d  increases  Psph  increases. 

.".  When  b/sin  9  >  d,  Psph  increases  as  d  increases. 

Since  when  b/sin  9  >  d,  Psph  T  as  d  T 

and  when  b/sin  9  <  d,  Psph  1  as  d  T  , 

the  Psph  must  reach  its  maximum  when  d  =  b/sin  9,  or  when 

d  =  [EMD  +  W/2]/sin  9. 

This  argument  is  readily  extended  to  sticks  containing  more  than  two 

weapons  by  arguing  that  each  interval  should  be  set  so  as  to 

produce  the  maximum  Psph,  hence  each  interval  between  weapons  should 

be  set  so  that 

INT  =  2d  =  2 [EMD  +  W/2]/sin  9 


and 


or 


SL  =  (n-1)  INT  =  2(n-l)  [EMD  +  W/2]/sin  9 

INT  =  [2  EMD  +  W]/sin  9 

SL  =  (n-1)  (2  EMD  +  W)/sin  9. 
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APPENDIX  B 

LIMITS    OF    INTEGRATION   RAILROAD/ROAD  TARGETS 
In   order    to  discuss    the    limits    of    integration   define    two   co- 
ordinate   systems,    a   X-Y  coordinate   axis    parallel   and    perpendicular 
to   the    target   dimensions    and   a    R-D   coordinate    system    (used    to  measure 
pilot   error)    parallel   and   perpendicular   to   the    flight    path.      See 
Figure    15. 


Coordinate   axes    in   railroad/road    problem 
FIGURE    15 
Once   again   in   order    to   reduce    the    notation  assume    that    n  =  2    (i.e. 
the   stick   consists    of   two  weapons).      Furthermore,    for   any  given 
approach   angle    (0),    set    the   bomb    interval    (INT)    for    these   weapons 
at    the    optimal   distance, 

INT  =    (2   EMD  +  RRW)/sin   9. 
Now,    there    is    a    target  kill    if   the   CS    falls   within  certain    limits. 
See   Figure    16. 
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Limits  of  position  of  stick  center 
needed  to  produce  a  target  kill 

FIGURE  16 

where  y  =  upper  limit  for  CS  and 

y  =  lower  limit  for  CS . 

The  question  then  is  what  is  the  mathematical  representation  of  these 

limits. 

From  Figure  16  "it  is  apparent  that  y  =  b  +  d  sin  0  and  that 

y.    =  -y  .   From  Appendix  A 


b  =  (EMD  +  RRW/2) 

d  =  (EMD  +  RRW/2 )/s in  9 


hence 


y  =  (EMD  +  RRW/2)  +  (EMD  +  RRW/2)  sin  9/sin  0 
=  2  EMD  +  RRW 


and 


y      =   -    [2    EMD  +  RRW]; 
however,    the    pilot   miss-distance    is   measured    in   the   R-D  coordinate 
system  and   must   be    rotated    into   the   X-Y  coordinate   system  before    the 
integration  can   take    place.      Therefore,    introduce   a    change    of  variable 
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R  =  X  cos  9  +  Y  sin  9 
D  =  -X  sin  9  +  Y  cos  9, 
The  Jacobian  of  this  transformation  is 


d(R,D)|         ("cos    9      sin   9]   =    . 
cKX,Y)|      "  L-sin   9   cos    9J 


and    the 


Psph 


a,      [2   EMD+RRW] 

J    r 


J    J 


fx  y(X,Y)    dXdY 


[2   EMD+RRW 


where 


fx      (X,Y)    = 


n T~   exP 

2TTaR°D 


1  — 


■  I  s  a  S  ] 


and 


A  = 


2       2 

cos  9  sin  9  -cos9  sin9  ,  cos9  sin9 

2       2'       2  2 

°R      CTD        QR  °D 


2       2 
cos9  sin9   cos9  sin9   sin  9   cos  9 

2  2        2       2' 


UR  UD       UR       R   /• 

The  double  integral,  although  having  limits  of  integration  that 
are  easy  to  evaluate,  has  an  integrand  that  precludes  a  closed  form 
evaluation.   However,  this  integral  does  provide  some  insight  into 
the  problem.   Notice  that  the  limits  of  integration  are  independent 
of  the  approach  angle  (9).   This,  of  course,  means  that  no  matter 
what  value  9  takes  on,  the  infinite  strip  over  which  the  integration 
will  take  place  will  be  the  same.   Furthermore,  if  this  infinite  strip 
is  represented  by  coordinates  in  the  R-D  coordinate  system,  it  will 
remain  unchanged  as  far  as  size  and  shape  is  concerned.   These  facts 
are  true  because  the  transformation  involved  is  simply  an  area 
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preserving  rotation  of  axis.   This  can  be  seen  by  evaluating  the 
Jacob ian  of  the  transformation 

X  =  R  cos  0  +  D  sin  9 

Y  =  -R  sin  0  +  D  cos  9 


ld(X,Y)l      .  f"   cos   9  sin  §"]  ,  . 
lcUR,D)|     "  L-sin  9  cos   9J 


and   recalling   that    if   the   absolute   value   of   the   Jacobian   of  a    trans- 
formation equals    1,    then   the    transformation  does    not    change    the    scale 

4 
of  any  dimensions. 

The    problem  now  then   is    to   find    for  a   given  bi-variate    normal 

distribution   the   orientation   of  a    fixed  dimension   infinite   strip 

so  as    to   sweep   out    the   maximum   probability.      See  Figure    17. 


Contours    of   equal   probability   of  a 


bi-variate    normal  distribution,    an  >  O 

R  D 

FIGURE    17 

Obviously,    the   maximum  probability   is   swept    out   if   the    infinite    strip 

lies   along    the   major  axis    of   the   ellispes    of  equal   probability  with 

its    center    line    coinciding  with    the   R-axis .      See   Figure    18. 


4 
Kaplan,   Wilfred,   Advanced  Calculus,    p.    201,   Addison  Wesley,    1957 
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7>-R 


Infinite    strip   sweeping    out   maximum   probability 
for   a   given  bi-variate    normal   distribution 

FIGURE    18 

However,  when  9=0   is  used  to  evaluate  the  bomb  interval  where 

INT  =  [2  EMD  +  RRW]/sin  0 

the  interval  >  ». 

The  meaning  of  this  result  is  that  the  interval  should  be  infinitely 

long.   This,  of  course,  is  impossible  to  achieve  but  it  does  point  out 

that  the  bomb  interval  should  be  as  great  as  possible,  while  still 

satisfying  the  conditions  for  interval  optimality.   Recall  that  there 

has  been  defined  a  maximum  and  minimum  release  altitudes.   These  two 

constraints  limit  the  bomb  interval.   See  Figure  19. 

Clearly  SL    ~  [ALT    -  ALT  .  I  cotan  a  where  a  =  dive  angle, 
max       max      min 


and  since 


INT  =  SL/(n-l) 


INT 


max 


(ALT    -  ALT  .  )  cotan  a 
max min 

(n-1) 
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MIN  RELEASE  ALT 


MAX  RELEASE  ALT 


approx, 


Maximum  stick  length,  n  =  4 

FIGURE  19 

When   n  =  2    this   expression   reduces    to 

INT  =    (ALT  -  ALT    .    )   cotan  a. 

max  max  mm 

Now   the    interval  must   also   satisfy   the   expression 

INT  =    [2   EMD  +  RRW]/sin   9 

max 

in   order    to   produce   an   optimal   stick    length.      Hence    it   is    possible 
to   determine    the   best    possible   approach  angle: 


(ALT  -  ALT    .    )   cotan  a 


max  min 


n-1 


=  INT  =    (2   EMD  +  RRW)/sin  0 

max 


or 


sin   9   = 


(ALT          -  ALT    .    )    cotan  a 
max min 

(n-1)    (2   EMD  +  RRW) 


or 


opt 


=  sin         r 


1    j[ALT 


-  ALT    .     ]    cotan  a 
max  min 


(n-1)    (2   EMD  +  RRW) 
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In   order    to    provide    specific    results   and    to   satisfy    the    requirements 

of    the   computer    programs    the   value    for      (ALT  -   ALT    .     )    cotan  ct| 

K  K     °  L  max  mm  -I 

was    arbitrarily   set    equal    to    1000    feet.      This    corresponds    to  a    maximum 
and   minimum   release   altitude    difference    of  approximately    1000'    if    the 
dive   angle   =  45    . 

Although    it    is    now   possible    to   determine  what    the    best   approach 
angle    is,    it    is    not    possible    to   determine  what    Psph    is    for   any 
specific   value    of   9   nor    is    it    possible    to   evaluate   how  much    the    Psph 
is   degraded    by   varying    from   the    optimal   conditions.      To  achieve    these 
results    it    is    necessary    to   evaluate    Psph    for  all    9. 

Since    rotating    the    pilot   miss-distance   distribution  coordinate 
system   into    the    target    coordinate    system   produced   a    difficult    inte- 
grand,   the    reverse    procedure   might    produce    better   results.      That    is, 
rotate    the    limits    of   integration    into    the   R-D   coordinate    system. 
Recall   that    the    integrals    involved   are 


°°     2EMD+RRW 
Psph   =  |     j 

-oo   -[2EMD+RRW; 


fx  y(X,Y)dYdX 


and 


Psph   - 


/( 


fR  D(R,D)dRdD 


appropriate 
strip 


where 


fX  Y    (X' 


Y)=2^e*p[-I^A    (y) 


and 


fD   n    (R,D)    =  0_  -     exp 


R  D 


2mRaD 


/  2  T 

i/r         D 


2       2"        2J 
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consider  the  transformation 

X  =  R  cos  9  +  D  sin  9 

Y  =  R  sin  9  +  D  cos  9 

once  again 

d(X,Y)   = 
d(R,D) 

The  problem  then  is  to  use  this  transformation  on  the  limits  of 

integration.   Recall  what  the  area  of  integration  under  consideration 

looked  like  in  the  X-Y  coordinate  system.   See  Figure  20. 

P3  /K  P4 


* 


7^ 


A  _ 


Area  of  integration  in  X-Y  coordinate  system 
FIGURE  20 
Since  the  transformation  used  is  area  preserving  all  that  is 
necessary  is  to  represent  this  infinite  strip  in  the  R-D  coordinate 
system. 


Define  points 


PL  =  (X1}  Yx) 
P9  —    (X~  ,  Y9 ) 


P«  —  (Xq ,  Y„) 
p4  "  <X4>  Y4> 
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Then  use  the  transformation 


Pl  (X1'  V  ~*   Pl  (R1'  V 
P2  (X2,  Y2)  -  P2  (R2,  D2) 


P4    <V    V   -   P4    (R4'    V' 
Now   since    (P~    -    P1 )   defined   a    line    in   the   X-Y  coordinate    system,    it 
defines   a    line    in    the    R-D  coordinate    system  and    the   equation    for    this 
line    has    the    form: 

DR  =  mLR  +  bx 

where    the    subscript   R  refers    to    the    right    hand    line    in  Figure    21, 

and 

R2    "  Rl 

m..    =   slope    = 


1        —        D2    -    Dl 


R.    -R, 


'1   =   ^tercept   =  R±    -  Q    _   J)  ^ 


2      "r 

Likewise    the    points    (P,    -    P   )    define   a    similar    line 

D      =  m2R  +  b 


where 


RA   -R3 


R,    -   R 


G^: 


b2   '  R3       VJ).    -   X>J  "3 


These   equations    produce   an  area    of    integration   of    the    form   seen    in 
Figure    21. 


Sisam,    C.    H. ,    College  Mathematics,    p.    185,    Holt,    1957 
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Area    of   integration   in  R-D  coordinate    system 
FIGURE   21 


Then 


00     D      =  m2R+b2 


Psph   = 


o exP 

2ttgrgd 


[-*(* 


27  J 


dDdR, 


-°=  D      =  m  R+b 
R  11 


This    integral   is   easily  evaluated   on  a   computer   by  a    numerical 
integration   technique    (see   Appendix  D)   allowing  R  to   range    from 
-oo    to  +oo. 

It    is    now  possible    to   evaluate    Psph  at  any  approach  angle    (9), 
determine    the    optimal  approach  angle,   evaluate    the    optimal   bomb 
interval   and    furthermore,    determine    the   effects    on  Psph   of  changing 
from   the    optimal    tactics. 

Although   this   development   was    performed  when   n  =  2    it    is   easily 
extended    to   n   =   3,    4,    ...,    n.      The    only    portions    of    the    development 
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that  change  are  the  formulae  for  determining  y   and  v..   Recall  in 
o  '  u       I 

the  case  n  =  2 

y   =  b  +  d  sin  9 
y  u 

y  =  -[b  +  d  sin  9] 
and 

b  =  EMD  +  RRW/2 

d  =  [EMD  +  RRW/2 ]/sin  9. 
For  a  stick  of  greater  than  n  =  2, 

y   =  b  +  (n-1)  d  sin  9 

yl  =  "yu 

This  is  most  easily  seen  by  determining  geometrically  just  where 

the  stick  center  must  be  positioned  in  order  to  just  produce  a 

kill  (y.  is  depicted  simply  for  ease  of  construction).   As 

illustrated  in  Figure  22  there  is  a  kill  when 

y,  <  CS  <  y  . 
I  —    —  u 


The  d 
but 

and 


•4.      •    i   j  •    i    i      (n-1) 
lstance  involved  is  clearly  -  -y-*— - — <- 


1 


(INT)/sin  9 J  +  b 


J 


INT  =  2d 


=  -  -l(n-l)  d/sin  9 J  +  bl 

f  1 

^(n-1)  d/sin  9  +  b,  . 
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Lower  limit  of  CS  producing  a 
target  kill  n  =  6 

FIGURE  22 
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APPENDIX   C 

LIMITS    OF    INTEGRATION   BRIDGE   TARGETS 
For    ease    in    notation  assume    n   =   2,    and    9  >   (?  .      The    optimal    bomb 
interval    is    then  determined    in   exactly    the    same    manner   as    for    the 
railroad/road    problem.      That    is 

INT   =    (2    EMD  +  BW)/sin   9   =   EBW/sin    9. 
Next,    the    question   of  where    the    center    of    the    stick    can    fall   and 
still    produce    a    kill    is    answered.      Clearly    until    the    EBL  becomes    a 
factor,    y     and   y     are    the    same   as    for   the    railroad/road    problem, 


name  ly 


BW 
y      ■  ~  +  EMD  +  d   sin   9 
u        2 

EBW    ,    EBW  sin   9   _  _-„ 
—      ~       r      x  ;       —   -   E.BW 

2  2      sin    9 

y         =      _y         =      -EBW. 

J  I  ■'u 

However,    the   effective    bridge    length   affects    the   area    of   integration 
by   defining   right   and    left   hand    limit    points    in   the    X-direction. 
Consider   y    ,    the    left   and    right   hand    limits   are   determined   by   the 

most   extreme    position   of  CS  when   the    forward   most  weapons    just  kills 

R  R  R 

the    target.      See   Figure   23.      These    two   points,    namely   P..       (X,     ,    Y,    ) 

and    P.       (X1     ,    Y,    )   where    the   superscript  denotes    left    or    right    hand 

limit,    can  be    represented   by 


[S£  +  d   cos    9" 


["EBL        EBW  cos    9 
L    2  2      sin   9J 


["EBL        EBW 
L  2  2 


:otan   9 
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v  R   EBL    ,      _    EBL   EBW    ,.    _ 
X   =  — y-   -   d  cos  9  =  —r-   -   — ~—  cotan  9 


^X 


CI  of  bomb 


Right  and  left  limits  of  y. ,  n  =  2 

FIGURE   23 

R  L 

These    two    limits    in   the   X-direction,    X..    ,    X..      are    then  active   until 

the   CS   coincides   with   the    center    line    of   the    target.      The    limits    of 
integration   then  change.      See   Figure   24. 


^X 


Limits  of  integration,  CS  coinciding  with 
center  line  of  target 

FIGURE  24 
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Now 


— \-   d  cos  9  I 


"EBL  ,  EBW        _~| 
j  — ^—  +  -5—  cotan  9j 


Y„   =  0 


EBL 


+  d  cos  9 


EBL  J  EBW    .    Q 
— 1 —   cotan  9 


0 
EBL 


EBW 


cotan  9 


Y„   -  0 


EBL  ,  EBW    ^    _ 
— = 1 r—  cotan  9 


R     L 
The  limits  X„  ,  X„   hold  until  CS  falls  on  y  and  this  completes  the 
j    3  u 


polygon.   See  Figure  25 
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Area  of  integration  in  X-Y  coordinate  system 
FIGURE  25 
P4    consists  of  X^L  =  X?L,  y/  =  Yu  =  EBW;  and  P^  consists  of 


4     3  '  x4 
v  R   v  R  _  R 
X4  =  X3  '  Y4  =  yu  =  EBW 


R 


The  polygon  is  now  well  defined  by  the  eight  points  (P  L, 
,  ...,  P^  )  and  by  using  the  transformation 

R  =  X  cos  9  +  Y  sin  9 
D  =*  X  sin  9  +  Y  cos  9 
the  same  polygon  can  be  represented  in  the  R-D  coordinate  system. 
See  Figure  26. 

The  problem  is  now  essentially  the  same  as  the  railroad  problem 
That  is  to  say 


.'41. 


Psph  = 


// 


area  of 
polygon 


o exp 

2mRGD 


2  2 

i  (£-  +  si: 

2  V     2  T        2. 
GR  °D 
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and  can  be  easily  integrated  by  numerical  techniques.   See  Appendix  D 
However,  the  limits  of  integration  need  to  be  systematically  defined. 


Area  of  integration  in  R-D  coordinate  system 

FIGURE  26 
The  procedure  used  is  to  define  the  right  and  left  hand  limits  of 
integration  for  the  deflection  error  (D)  by 

DT  =  m.R  +  b.    j  =  1,  ....  4 

L    J      J 

D„  =  m.R  +  b.    i=l,  . . . ,  4 
R     L       1 

where  m^  and  b  are  the  appropriate  slope  and  intercept  for  the 
portion  of  the  boundary  of  the  polygon  involved  in  the  integration, 
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The   values    of  m„  and   bv  can  be   calculated    from   the    same    formulae   as    in 
K  K 

the    railroad/road    problem. 
Therefore  , 


DT=m.R+b. 
00      L     J         J 


Psph 


-l\ 


-°°  D   =m.R+b. 
R      1         1 


2mRaD 


exp 


[-K; 


2  2 

2  2J  J 

GR  GD 


dDdR 


i   =    1,    ....    4 
J    =    1,    -..,    4; 

and    the   appropriate   m     and   b     are    used    for    the   value    of    dR  during    the 

K       K 

numerical  integration. 

One  interesting  note  is  that  when  9  =  j?  the  polygon  collapses 
somewhat  since 


cotan  9  =  cotan  ?  = 


EBL 
EBW 


Hence 


EBL   EBW  EBL 
.  2      2   EBWJ 


=  -  EBL 


R   EBL   EBW  EBL 
1     2   "   2   EBW  " 


X„   -  -  EBL 


X„   =  0 


X, 


EBL   EBW  EBL 
2      2   EBW 


=  0 


See  Figure  27. 


R   EBL   EBW  EBL  _ 
X3   "  2  +  2   EBW  "  EBL 


X.L  =  0 
4 


X.   =  EBL. 
4 
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Area  of  integration  X-Y  coordinate  system, 
approach  angle  equals  critical  angle 

FIGURE  27 

Now  suppose  that  the  approach  angle  is  less  than  p,  this  type  of 

attack  is  exactly  the  same  except  the  roles  of  EBL  and  EBW  are 

reversed.   The  optimal  bomb  interval  is  now 


INT  = 


EBL 


cos  0 

and  the  formulae  for  the  8  points  in  the  polygon  are  exactly  the  same 

except  EBL  is  replaced  by  EBW  and  vice  versa  and  due  to  the  manner  in 

which  9  is  measured  cotan  9  is  replaced  by  tan  9,  see  Figure  28. 

Y 

t 


Area  of  integration,  9  <  |3 
FIGURE  28 


65 


Once  the  polygon  has  been  described  the  remaining  development  is 
exactly  the  same  as  when  9  >  3  and,  therefore,  is  not  discussed. 
Notice,  however,  that  when  9=3  the  formulae  reduce  to  those  obtained 
previously  when  9=3  (as  they  should).   One  implicit  assumption  in 
this  development  is  that  the  target  is  assumed  to  be  only  of  length 
BL.   That  is  to  say  cratering  the  approaches  to  the  bridge  is  con- 
sidered a  miss.   This  type  of  damage  can  easily  be  accounted  for 
however  by  simply  changing  BL,  or  changing  the  lethality  function. 

The  next  problem  is  to  extend  this  development  to  an  arbitrary 
number  of  weapons.   Suppose  9  is  greater  than  3,  note  that  for  each 
weapon  in  the  stick  there  are  4  limit  points,  directly  analogous  to 
the  four  corners  of  the  target.   The  polygon  over  which  the  integration 
is  performed  then  must  have  4n  corners  for  any  n.   Note  that  the  opti- 
mal bomb  intervals  are  equal  for  a  given  9.   That  is 

INT  =  EBW/sin  9 

and  y.  and  y  are  defined  in  the  same  manner  as  in  the  railroad/road 
I      u 

problem 


where 


hence 


and 


y   =  b  +  (n-1)  d  sin  9 

u 


EBW        EBW 
b  =  — r- ,  d  = 


2  '     2  sin  9 


EBW  ,  .      _.  EBW  sin  9 
u    2  2   sin  9 


y,  =  -y 


u 
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Now  the  right  and  left  hand  limit  points  for  y  must  be  determined 
They  are  defined  by  the  most  extreme  position  of  CS  when  the  most  for- 
ward bomb  just  produces  a  kill.   See  Figure  29. 


Limits  of  integration  for  y  ,  n  =  6 
FIGURE  29 
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Clearly 


L  =    _    |~EBL 


V      =r +  ((!  -  0 


INT  +  d 


l)  cos    Oj 


,"> 


--[™  +  ((ai^)[2d]  +  dj 


] 


j-SL 


+   (n-1)   d   cos    9 


and 


XiR  =  Mi   -     (n-l)    d   cos    0. 


1  2 

In  a  similar  manner  each  point  of  the  left  and  right  hand  limits  can 
be  constructed,  using  the  following  formulae: 

„  L    „  L     JEBL  ,  ,    N  EBW   t    J      _i   o   q       0n    i 
X.  ,  .  =  X.   =  -  <-=-  +  (n-s)  -57-  cotan  9  \     ,   s  =  l,  3,  5,  .  .  .  ,  Zn-1 

1+1      1  I  *■ 

i=l,  3,  5,  . . . ,  2n-l 
YXL  =  \   EBW 


Y'J.1  =  Y-    =  9  EBW'  P  =  -(n-2)'  "(n-4)» 

1=2,  3,  . . . ,  2n-l 

Y9L  =  ?  EBW 
2n    2 


,  -2,  0,  2,  ...  n-2. 


X. ,.  =  X.R  = 

l+l    i 


EBL    .        N  EBW 
~2 (n-s)  -j-  c 


otan  9  , 


s=l,  3 ,  . . . ,  2n-l 
i=l,  3,  . . . ,  2n-l 


Y.R  =  Y.L,   i=l,  2,  3,  ...,  2n. 


See  Figure  30. 

The  formulae  when  9  <  |3  once  again  are  identical  except  for  the 
replacement  of  EBW  by  EBL  and  vice  versa  and  cotan  9  by  tan  9.   Once 
the  polygon  for  an  arbitrary  n  has  been  defined  in  the  X-Y  coordinate 
system  the  remaining  development  is  identical  to  the  case  for  9  >  0. 
The  only  problem  arising  being  the  "bookkeeping"  of  the  proper  left  and 
right  hand  limits  in  the  numerical  integration. 
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Area  of  integration  for  n  =  6 
FIGURE  30 
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APPENDIX  D 

NUMERICAL  EVALUATION  OF  THE  PROBABILITY  INTEGRAL 

In  order  to  discuss  the  numerical  integration  method  for  a  given 

set  of  parameter  values  (i.e.,  target  dimensions,  o  ,  a  ,  9,  n,  EMD, 

R    D 

g  (9))  the  optimal  stick  length  and  bomb  interval  is  determined.   The 
y 

area  of  integration  in  the  X-Y  coordinate  system  is  then  defined 
using  the  formulae  developed  previously.   This  polygon  is  then 
rotated  into  the  R-D  coordinate  system  by  the  appropriate  transform- 
ation and  the  following  type  of  numerical  integration  is  performed. 
Recall  the  integral  to  be  approximated 

Psph  =    if    ^V  exp  [-  £  (rr  +  S)]- 

-         D  R  a„  CL. 

area   of  R  D 

polygon 
This  integration  is  approximated  by 


D    /•*  RT  °L 

8a„/3  r\   ri 


Psph  = 


1=1  J  b  Jr     R  °  °R     °D 


R.  D 

8aR/3 

where   2  represents  the  addition  of  small  strips  of  probability  and 

c=l 

8aD/3      determines  the  number  of  strips.   For  example,  suppose  the 
R 

area  of  integration  in  the  R-D  coordinate  system  is  as  seen  in  Figure 
31. 

First,  the  lower  limit  R  of  integration  is  determined  and  R  .  is 

J-i  J. 


set    equal    to  R 

Li 
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Area  of  integration  in  R-D  coordinate  system 
FIGURE  31 
The  b  denotes  lower  value  of  the  AR..  interval  and  1  denotes  the  first 
interval  in  the  summation.   The  upper  limit  of  AR   is  then  defined  by 

R^  =  Rb  +  3. 
Next  the  center  point  R*  of  AR   is  determined  by 

\  +  hh 

and   the    right   and    left   hand   end   points    of   D  are   determined: 

D,R  -  m.R*  +  b. 
Ill 

D,L  =  m.R*  +  b., 
1  J  J 

using  the  proper  equation  for  the  boundary  of  the  polygon. 
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1_  «T« 

Next    R        and   R        are    normalized   and    the 

-R,T  Rb 


<  R  <  is   computed   by 


PROB 

°R  WR 

RlT/aR  2 

IK  u 

r     i    -  — 

PROB    (1)    =        J  2^  e  du. 

Rl   /aR 

R  L 

Next   D,      and   D,      are    normalized   and 

•  D,        D 


PROB  —  <  D  <  

L0D  "   -°D 


is  computed  by 


D1L/0D     _  J 
PROB  (2)  =   J    |^  e  2  du. 


D1R/°D 


Next  [PROB  (2)]  x  [PROB  (1)]  is  computed  and  saved. 

The  program  then  sets  R„   equal  to  R   and  repeats  the  process 
adding  the  increase  in  probability  to  the  previous  total.   This  pro- 
cedure is  repeated  until  the  uppermost  limit  of  R  is  reached  at  which 
time  the  total  value  of  the  probability  accumulated  is  printed.   The 
entire  procedure  is  repeated  for  each  0. 

One  further  comment,  in  the  railroad/road  problem,  since  R 
theoretically  varies  from  -«  to  +°°,  the  integration  is  limited  to 
±4oR. 

The  procedure  is  exactly  the  same  when  EAA  is  present  but  the  area 
in  the  X-Y  plane  is  slightly  different.   No  problems  are  encountered 
when  9  <  9  .   That  is  as  long  as  the  bomb  interval  used  is  less  than 
or  equal  to  the  optimal  bomb  interval  for  the  9  in  the  integration; 
however,  when  the  bomb  interval  used  exceeds  the  optimum  bomb  interval 
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(i.e.  9  >  0n)  areas  appear  where  the  CS  may  fall  and  no  kill  occur, 
Figure  32  illustrates  this  condition  for  n  =  6. 


Area  of  integration  EAA  present,  9  >  0 
FIGURE  32 
Formulae  for  the  construction  of  this  area  are  easily  determined  but 
since  the  computer  program  developed  is  a  special  case,  that  is  the 
number  of  bombs  is  constant  at  6,  general  formulae  are  not  included. 
The  problem  was  investigated  only  because  questions  arose  about  the 
validity  of  using  errorless  approach  angles  but  since  the  results 
obtained  were  not  greatly  different  than  when  9  was  assumed  perfectly 
known  the  investigation  was  discontinued. 
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APPENDIX  E 


AIR  WING 


GLOSSARY 
Tactical  unit  aboard  aircraft  carrier  made  up  of 
approximately  6  squadrons.   Normally,  the  largest 
unit  of  aircraft  that  operate  together. 


APPROACH  ANGLE 


Acute  angle  measured  from  center  line  of  target 
to  horizontal  projection  of  flight  path  during 
dive . 


BOMB  INTERVAL: 


Horizontal  distance  measured  along  flight  path 
between  the  centers  of  impact  of  two  adjacent 
bombs . 


BOMB  STICK: 


Line  of  bombs  dropped  at  programmed  intervals  in 
a  single  dive.   Normally  falling  in  a  straight 
line . 


DIVE  ANGLE 


Acute  angle  measured  in  vertical  plane  between 
the  horizontal  and  the  flight  path  during  the 
dive  . 


EFFECTIVE 
MISS -DISTANCE 


Maximum  distance  any  target  dimension  may  be  from 
the  center  of  impact  of  the  weapon  and  still  pro- 
duce a  kill. 


FLIGHT 


Several  airplanes  operating  together  as  a 
tactical  unit. 
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GENERAL  PURPOSE      Conventional  ordnance  constructed  of  high 
BOMBS:  explosives  and  iron  alloys.   Produced  in  different 

weights  and  used  on  ordinary  targets. 

KILL  THE  TARGET:     With  respect  to  a  railroad/road  target,  kill 

implies  cratering  the  road  so  as  to  make  it 
temporarily  impassable;  for  a  bridge,  collapsing 
the  bridge  is  considered  a  kill. 

LIGHT  JET  ATTACK     Turbo-jet  aircraft  designed  primarily  to  carry 
AIRCRAFT:  air  to  ground  weapons.   Much  smaller  in  size  than 

strategic  bombers  such  as  B-52 ,  B-47 . 

PILOT  MISS-DISTANCE   Probability  distribution  describing  pilot  induced 
DISTRIBUTION:        errors  between  the  intended  point  of  impact  of  a 

weapon  and  its  actual  point  of  impact. 

ROLL-IN:  Initiation  of  the  tracking  dive  in  a  bomb  run. 

SINGLE  PASS  HIT:     Killing  the  target  in  one  pass. 

SQUADRON:  Tactical  unit  aboard  aircraft  carrier  made  up  of 

approximately  15  aircraft  of  the  same  type. 
Smallest  organizational  unit. 

STICK  BOMBING:       Tactic  involved  in  dropping  several  bombs  during 

a  single  dive  using  programmed  intervals. 

STICK  LENGTH:        Horizontal  distance,  measured  along  horizontal 

projection  of  flight  path,  from  center  of  impact 
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of  the  first  weapon  released  to  center  of  impact 
of  last  weapon  released  during  a  single  run. 

TRACKING  DIVE:       Dive  in  which  bombs  are  intended  to  be  dropped. 

Pilot  attempts  to  correct  for  all  known  errors 
and  to  place  the  geometric  center  of  the  bomb 
stick  on  the  geometric  center  of  the  target. 
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